Quantum degeneracy and interaction effects in spin-polarized Fermi-Bose mixtures 
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Various features of spin-polarized Fermi gases confined in harmonic traps are discussed, taking 
into account possible perspectives of experimental measurements. The mechanism of the expansion 
of the gas is explicitly investigated and compared with the one of an interacting Bose gas. The 
role of interactions on the equilibrium and non equilibrium behaviour of the fermionic component 
in Fermi-Bose mixtures is discussed. Special emphasis is given to the case of potassium isotopes 
mixtures. 

I. INTRODUCTION 



After the achievement of Bose-Einstein condensation 
(BEC) in dilute vapors of alkali atoms confined in mag- 
netic traps Jl|-§|, a challenge for future experiments is 
the cooling of samples of fermionic atoms, down to the 
degenerate regime. 

This perspective has already motivated several theoret- 
ical works (see, for example, [||J5)) devoted to the study 
of a possible superfiuid phase occurring at very low tem- 
perature, analog to the one exhibited by superfiuid 3 He 
and superconductors . The behaviour of this transition 
depends crucially on the sign and the size of the scat- 
tering length and various predictions have already been 
made to estimate the value of the corresponding critical 
temperature. Another peculiar phenomenon exhibited by 
finite Fermi gases at low temperature is the occurrence 
of quantum shell effects [||, bringing interesting analo- 
gies with the physics of atomic nuclei and metal clusters. 
The observation of such phenomena is seriously limited 
by the difficulty of reaching the required regime of very 
low temperatures. In fact both superfluidity and shell ef- 
fects are expected to occur at temperatures much smaller 
than the Fermi temperature. For a realistic perspective 
of experimental measurements in the immediate future 
it is consequently useful to explore quantum degeneracy 
effects taking place at higher temperature. These effects 
have been already the subject of several theoretical works 

HBJ§. 

A crucial role in the experimental achievement of the 
low temperature regime needed to observe quantum de- 
generacy effects is played by evaporative cooling. This 
technique requires frequent collisions in order to ensure 
fast thermalization. Scattering between spin-polarized 
fcrmions in s wave is inhibited by the Pauli exclusion 
principle. As a consequence one should look for alter- 
native scattering processes involving different spin states 
or atoms belonging to different species. Fermi-Bose or 
Fermi-Fermi mixtures are possible candidates. The re- 
sulting scattering processes will eventually provide the 
relevant thermalization mechanism, allowing for the cool- 



ing of the sample (sympathetic cooling |15|,|16|] ) . 

In the following we will mainly discuss the case of 
Fermi-Bose mixtures where the effects of two-body inter- 
actions can be significant also on the equilibrium proper- 
ties of the fermionic component because of the occurrence 
of Bose-Einstein condensation which strongly enhances 
the value of the central density below the critical temper- 
ature, thereby emphasizing the effects of the mean field 
interaction. In this context we will limit the discussion 
to the case of spin-polarized Fermi gases where the inter- 
action between fermions can be safely neglected. Notice 
that in the case of unpolarized Fermi gases also the in- 
teraction between fermions can produce significant effects 
on the equilibrium properties if the scattering len ght is 
large as happens, for example, in the case of 6 Li PJlq]. 

The aim of the present paper is to discuss some of the 
key features exhibited by spin-polarized trapped Fermi 
gases, with special emphasis to the dynamics of the ex- 
pansion following the switching off of the trap (Sect. ||) 
and the role of interactions in the case of Fermi-Bose mix- 
tures (Sect. pi) ). We apply the results to the case of 40 K, 
the fermionic isotope of potassium, which has been re- 
cently cooled in a magneto-optical trap JT^ j . Potassium 
is in fact a good candidate for the experiments proposed 
here because of the presence of two bosonic isotopes, 39 K 
and 41 K in addition to the fermionic one. The possibility 
of magnetically trapping these atoms has been already 
demonstrated flS|] . Collisional cross sections for potas- 
sium have been calculated in p^,pO| . In the following 
we assume the values given in |}20|. The analysis pre- 
sented in the present paper can be easily extended to 
other fermionic atoms which are presently investigated 
experimentally as, for instance, 6 Li pifl. 



II. THE IDEAL FERMI GAS 

The thermodynamic behaviour of ideal trapped Fermi 
gas has already been investigated [^- 10 (| . In this section 
we summarize the main results and we derive explicit 
formulae for the expansion of the fermionic cloud after 
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turning off the confining potential. These results could be 
relevant in view of future measurements on the expanding 
cloud. 

Let us consider N fermions trapped in an axially sym- 
metric harmonic potential 



V{r±,z) = im (wiri + uj 2 z z 2 ) 



(1) 



and ignore the effects of the two-body interatomic force. 
The particle distribution function can be written, in semi- 
classical approximation, as 



n(r,k, (3) 



1 



1 



(2 7 r)3 exp(/ 3 ( ^ + ^( r± ^)- M/ )) + l' 



(2) 



where j3 = (KbT) 1 and fif is the chemical potential, 
fixed by the normalization condition 



N = / n(r, k, j3) dr dk 
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E 2 dE 



2{%uj) 3 J e^ E ~W + 1' 

(3) 

In eq. (||) uj = (cj^cl^) 1 / 3 is the geometrical average of 
the trapping frequencies. 

If more spin states are occupied, then (||) characterizes 
the distribution function of each component separately. 
It can differ from a component to another because the 
number of atoms in each spin state and the corresponding 
trapping potential are in general different. 

At T = eq. (||) allows one to calculate the Fermi 
energy defined by Ep = /U/(T = 0). One finds 



E F = K B T F = {6N) 1/3 huj. 



(4) 



The Fermi energy (Q) can be used to define typical length 
and momentum scales pp| characterizing the Fermi dis- 
tribution in coordinate and momentum space respec- 
tively: 



E F = -muiRi 



—vnw'Zi 



h 2 K F 
2m 



(5) 



where R± and Z are the radial and axial widths of the 
density distribution at zero temperature: 



"< r;T = 0) = Jnfe 



Rx 




while Kp is the width of the corresponding momentum 
distribution 



n(k;T = 0) 



8 N 
^K 3 



K 2 

J^p 



3/2 



(7) 



Eqs. @ and ([?]) hold for positive values of their argu- 
ments. Eq. (Q) is the analog of the most familiar mo- 
mentum distribution n(k;T = 0) = 3iV/(47r.K'f.)9(l - 
k 2 /Kp) characterizing uniform Fermi gases. 



At finite temperature the total energy E of the system 
can be written in the form 



E = 



1 



2(fno) 3 



E 3 dE 



(8) 



while the release energy, given by the energy of the sys- 
tem after switching off the trap is equal to E re i =15/2, 
because of the equipartition theorem. At low tempera- 
tures one has E ~ jNEp (l + ?f (^r)^ • In Fig- we 

show the release energy vs T for a gas of N fermions con- 
fined in a harmonic trap (solid line). In the same figure 
we also report the behaviour predicted by the classical 
gas (E re i — 3/2NKbT, dashed line) as well as the one 
of an ideal Bose gas with the same number of particles 
N, in the same confining trap (dot dashed line). For 
the ideal Bose gas the release energy is given by E re i — 
1.35N(K B T) 4 for T < T c where T c = OMUcjN 1 / 3 / K B is 
the critical temperature for Bose-Einstein condensation. 
Notice that, differently from the Fermi case, Bose gases 
exhibit a phase transition at T — T c . This temperature 
has the same N dependence as the Fermi temperature 
(see eq. (|^)) and is fixed by the same geometrical average 
of the trapping frequency. The general relation between 
the two temperatures is given by 



T F = 1.9T C 



N, 



1/3 



(9) 



where the suffix / (b) refers to fermions (bosons). The 
figure clearly shows that in order to observe effects of 
degeneracy in the release energy of a Fermi gas one should 
go to temperatures considerably smaller than the Fermi 
temperature. For example in order to obtain a deviation 
of 20 percent from the classical value one should work 
at T ~ 0.3 Tp. This represents a major difference with 
respect to the thermodynamic behavior of a Bose gas 
where the effects of Bose-Einstein condensation in the 
release energy show up immediately below the critical 
temperature as clearly shown in the same figure. 

Let us now consider the expansion of the cloud after 
turning off the trapping potential. To obtain the tem- 
poral evolution of the density profile one has to solve 
the Boltzmann transport equation. For cold and di- 
lute Fermi gases one can ignore two body collisions and 
the distribution function after switching off the trap 
will consequently follow the ballistic law n(r, k, f3, t) = 
no(r — — t, k, /3). where n is the distribution function 
at t = 0, given by (§). One can then easily calculate the 
time evolution of the spatial density during the expansion 
for which we find the simple analytic result 

n(r,(3,t) = 

6N f 1 1 

R\Z (ttEf/3)^ (1 + u 2 ± t 2 )(l + u 2 t 2 y/ 2 l3/2(Z) 1 ' 

where f n (z) = (Tin))- 1 J dy y n ~ x / \z~ l e y + 1) are the 
Fermi functions, S — exp(/3(^/ — V(r,t))) and 
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V(r,t) = -m 



1 + ojU 



2+2 



1 + uj 2 t 2 



plays the role of an effective potential fixing, at each in- 
stant, the shape of the distribution function. The chemi- 
cal potential has no time dependence, being fixed by the 
normalization condition ([?]). 

From ( |l0| ) one can extract the temporal evolution of the 
mean square radii of the system: 



>2 > = ^7^^r( 1 + w i <2 ) 

3N muj. 



3N mwi 



(11) 



(12) 



which have been written in terms of the release energy 
E re i = E/2 of the gas. The structure of these equations 
is independent of the temperature which enters the prob- 
lem only through the value of the release energy, fixing 
the initial value of the widths. Actually the release en- 
ergy fixes also the asymptotic behavior of the widths, as 
clearly shown by eqs. ( [ll"| - |T^ ). In Fig. |^ we compare the 
temporal evolution, at T — 0, of the root mean square 
radii of the Fermi gas (see eq. pT[]l^ ) with the one of 
a gas of interacting bosons. The Fermi gas corresponds 
to 10 6 40 K atoms initially confined by a harmonic trap 
with ujf = 100 Hz. The Bose gas instead corresponds 
to 10 6 39 K atoms interacting with a scattering lenght 
a bb = 80 ao 1 20 1 (ao is the Bohr radius) and initially con- 
fined by a magnetic trap with u) b = 100 Hz. The same 
parameters have been used to calculate the release energy 
of the interacting Bose gas as a function of temperature 
(dotted curve in Fig. [j]) in the framework of the theory 
developed in |2^| , Fig. [2] clearly shows that the Fermi 
gas expands more quickly than the Bose gas due to the 
significantly higher value of the release energy. 

Another useful quantity is the aspect ratio of the cloud 



R r (t) 




(13) 



Eq. (|13|) shows that for t — > 00 the cloud becomes spher- 
ical in shape even if it was initially strongly anisotropic. 
This behaviour is independent of the temperature and 
is due to the absence of collisions during the expansion. 
Notice that in a Bose gas the situation is quite differ- 
ent, the asymptotic distribution being anisotropic due to 
the occurence of Bose-Einstein condensation, as explic- 
itly pointed out in the first experiments on BEC |i]p| . 



III. FERMI-BOSE MIXTURES 

It is important to discuss how the scenario presented in 
the preceding section for the ideal Fermi gas is modified 
when the system interacts with a Bose gas confined in 



the same trap. As discussed in the introduction such 
mixed Fermi-Bose gases might become relevant in future 
experiments for the achievement of very low temperature 
regimes via sympathetic cooling. 

The study of Fermi-Bose mixtures has been already the 
subject of theoretical studies at zero JllJ as well as finite 
temperature |Q . Due to the occurrence of Bose-Einstein 
condensation the bosonic component is characterized by 
a high density in the central region of the trap where 
consequently interaction effects play an important role. 
The presence of the much more dilute Fermi component 
is not expected to influence the bosonic wave function in 
a significant way, so that the main effect of interactions 
between fermions and bosons will result in an additional 
external field acting on the Fermi component. The Fermi 
gas can then, in first approximation, be treated again as 
an ideal gas trapped by the effective potential 



V eff = V + g b fn b (r,T) 



(14) 



where V is the external potential (|l|) trapping the 
fcrmionic species, and the renormalization arises from 
the interaction with bosons. The parameter g b f — 
2irh 2 abf jm r is the Fermi-Bose interaction coupling con- 
stant fixed by the relative scattering length a b f and by 
the reduced mass m r = m b mf/(m b + rrif). The bosonic 
density n b can be calculated at thermal equilibrium us- 
ing standard procedures developed to describe Bose con- 
densed gases at finite temperature (see for example |22[ ) . 

The new equilibrium properties of the Fermi gas can 
be simply understood by looking at the shape of the 
potential ( |l4| ) as explicitly discussed in |l2|. A first 
important feature is that the interaction term in jl^ ) 
produces a weaker confinement. This results in an ex- 
pansion of the Fermi cloud with respect to the non in- 
teracting case and a consequent decrease of the aver- 
age density, as well as of quantum statistical effects. 
The effects are more pronounced at low temperature 
where all the bosons are in the condensate and one 
can use the Thomas-Fermi approximation to calculate 
the density distribution of the bosonic component. In 
this case one has n b = 15N b (Rl — r 2 )/ (8irRl) where 

R b = ahoimNhdhh/aho) 1 / 5 is the classical radius of the 
Bose gas, a bb is the boson-boson scattering lenght and 
a/10 = (ft/mbWf,) 1 / 2 is the oscillator length relative to the 
bosons. One finds that the effective potential felt by the 
fermions is given, for r < R b by ]12j 



eff 



V.. 2 



g bf m b uj b 
g bb rrifujj 



9bf_ 

.966' 



M6 



(15) 



where, for simplicity, we have considered an isotropic 
trap. For r > R b the effective potential instead coin- 
cides with the bare potential V. Eq. (|l^) shows that 
the oscillator constant is reduced by the interaction with 
bosons. 

As a first example we consider 10 4 40 K atoms and 10 6 
39 K atoms confined in the same harmonic trap (u>f = 
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u>b). This corresponds to a realistic case considering the 
present status of the experiments [|l7[|lg| l. For the inter- 
actions we make the assumption g^f ~ 6-5gbb pCp - 

The effect of the interactions is to push the fermionic 
cloud "outside" and to reduce the value of the central 
density. With the above choice of parameters we find, at 
T = 0, that the central density is reduced by ~ 35 per- 
cent, and the average square radius (r 2 ) increases from 
9.5 a\ to 11.5a 2 a . The importance of the interactions 
on the equilibrium properties of the Fermi component of 
the mixture is expected to depend in a crucial way on 
the temperature as well as on the relative concentrations 
of the two atomic species. For example at T — Tp one 
finds a much smaller effect due to the interactions. In 
this case the average square radius increases only by 2.5 
%, even if the central density is reduced by 25 %. Notice 
that in the configuration here considered the Fermi tem- 
perature is lower than the critical temperature for BEC: 
Tp = 0.4 T c . This shows that the effects of quantum de- 
generacy of the interactions can be observed only at very 
low temperatures. 

Concerning the role of the relative concentrations we 
point out that from an experimental point of view one 
can access quite a large range of concentrations. Indeed 
in JT^ ] 10 4 fermionic potassium atoms were magneto- 
optical trapped from a natural abundance sample. Since 
samples enriched by two orders of magnitude are avail- 
able, 10 6 40 K potassium atoms could be easily loaded. 
More recently more than 10 s bosonic isotopes have been 
trapped in a double MOT apparatus jl8| in the ultra high 
vacuum environmental condition required for evaporative 
cooling. This method can be easily extended to 40 K. 
Hence Bose-Fermi mixtures are expected to be available 
in the near future with a large variety of relative abun- 
dances. As a second example we have considered the case 
Nf = Nb = 10 6 . With these values the Fermi tempera- 
ture is larger than the BEC temperature (Tp — 1.9T C ). 
However the effect of the interactions are significantly 
reduced even at T = 0. For example the central den- 
sity has decreased only by 10 % with respect to the non 
interacting case. 

The above discussion suggests that the main features 
of the ideal Fermi gas discussed in section |l| are modified 
only in a semi quantitative way by the inclusion of the in- 
teractions with the bosons. Nevertheless the effects of the 
interaction could be observed experimentally by remov- 
ing the bosonic component from the trap at the end of 
the sympathetic cooling. This can be realized by optical 
or rf transitions, and corresponds to a sudden switching 
off of the interaction term at t = in the confining po- 
tential (|l5|) . The result is an undamped oscillation of the 
Fermi cloud, shown in Fig. ^, according to the law 



(r 2 ){t) = (r 2 ) + 6 cos(2u f t), 



(16) 



where (r 2 )o is the equilibrium mean radius of the 
fcrmions in the trap in absence of interactions, 5 is the dif- 
ference between the value of (r 2 ) calculated at t = and 



(r 2 )o- This oscillation, which is a direct consequence of 
the Fermi-Bose interactions, could be observed by imag- 
ing the atomic cloud during the expansion. 

Another consequence of the interaction between the 
two components is the occurrence of a shift in the fre- 
quency of the relative motion of the Fermi-Bose clouds, 
with respect to the predictions of the non interacting 
model. This can be easily calculated in the limit where 
the number of fermions is much smaller than the one 
of bosons. In this case the motion of the fermions is 
described by the potential 
of the fermionic cloud osci 

1/2 



and the center of mass 
lates with frequency uj = 



1 



9b f rn 



2 \ 



9bb rrif 

ployed above (u>f 



With the choice of parameters em- 

uj b , m b ~ rrif and g bf ~ 0.5g bb ) the 
frequency is reduced by about 30 percent. 

In conclusion we have shown that the behaviour of 
a spin-polarized degenerate Fermi gas, confined in har- 
monic traps, exhibits several interesting features, even 
in the absence of superfluidity. These features, which are 
the consequence of Fermi pressure effects, are directly ob- 
servable by looking at the expansion of the gas and are 
also sensitive to the interactions between fermions and 
bosons in atomic mixtures. 
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FIGURE CAPTIONS 

Fig. 1 Release energy in units of NKbTf as a function 
of temperature (in units of the Fermi temperature Tp). 
The solid line corresponds to an ideal Fermi gas trapped 
in a harmonic potential. The dashed line is the linear law 
of a classical gas, while the dot dashed and the dotted 
lines are the release energy for a Bose gas with the same 
number of atoms confined in the same trap in the non- 
interacting and in the interacting case respectively (see 
the text). 

Fig. 2 Temporal evolution in msec of the root mean 
square radii (in units of ato) of an ideal Fermi gas of 40 K 
(solid line) and of an interacting Bose gas 39 K (dashed 
line) at zero temperature after switching off the trap. 
The two curves refer to a trap with tu = 100 Hz and same 
number of particles N = 10 6 . The scattering lenght for 
the 39 K was taken from Q. 

Fig. 3 Oscillations at zero temperature of the mean 
square radius of 10 4 40 K atoms after the bosons (10 6 39 K 
atoms) are removed. Time and lenghts are in units of 
uj/tt and a^ Q respectively. 
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